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Abstract 

Background: Tubulin protein being the fundamental unit of microtubules is 
actively involved in cell division thus mabing them a potential anti-cancer 
drug target. In spite of many reported drugs against tubulin, fevw of them 
have started developing resistance in human p-tubulin due to amino acid 
substitutions. 

Methods: In this study we generated three mutants (F270V, A364T and 
0292E) using ModellerQvlO which were targeted with compounds from higher 
and lower plants along with marine isolates using iGEMD0CK2.0 to identify 
their residual interactions. 

Results: The mutant F270V does not bring in any increase in the binding affin- 
ity in comparison with the taxol-wild type due to their conservative substitu- 
tions. However, it increases the volume of the active site. A364T mutant brings 
a better binding among few of the marine and higher plants isolates due to 
the substitution of the non-reactive methyl group with the polar residue. But 
this leads to reduced active site volume. Finally the mutant 0292E from 
epothilone binding site brings a remarbable change in drug binding in the 
mutants in comparison with the wild type due to the substitution of un- 
charged residue with the charged one. But as such there was no change in the 
volume of the active site observed in them. 

Conclusion: Lower plants extracts were reported to exhibit better interactions 
with the taxol and epothilone binding sites. Whereas marine and higher 
plants isolates shows significant interactions only in the wild type instead of 
the mutants. In addition to this, the residual substitutions were also found to 
alter the conformations of the active sites in mutants 
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Introduction 

Tubulin being the fundamental unit of mi- 
crotubules is critically involved in chromo- 
somal segregation, cell division, motility and 
intracellular transportation \ They are made 
up of a and P subunits alternatively arranged 
in a lateral and longitudinal manner. The lat- 
eral contacts involve the interactions of Hl- 
S2 loop and helix H3 with the M-loop of the 



adjacent protofilament. Thus, 13 protofila- 
ments associate laterally and are found to be 
more electrostatic and less hydrophobic than 
the longitudinal contacts ^'^ . The alternative 
arrangement of a and P subunits results in 
longitudinal interactions which are classified 
into intra- and interdimer interfaces. The 
intradimeric interface is observed between p 
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and a subunits whereas the interdimericinter- 
face is found between a and (3 subunits 
Further, the longitudinal contact involves the 
interaction of H8 of a tubulin with H11-H12, 
T5, T3 and y-phosphate of the adjacent subu- 
nit. Similarly, T7 of a-tubulin shows interac- 
tions with phosphates T2, Tl, H7 and Gua- 
nine 

(3-tubulin subunit comprises three distinct 
domains which include the N-terminal do- 
main, intermediate domain and C-terminal 
domain. The N-terminal domain harbors the 
nucleotide and has 6 parallel beta strands (Sl- 
S6) along with the same number of alpha he- 
lix (H1-H6). The intermediate domain has 
strands S7-S10 with three helices H8-H10 
which accommodates taxol. The C-terminal is 
made up of two antiparallel helices H11-H12 
which interacts with Microtubule Associated 
Protein (MAP) 

Even though the overall domain architec- 
ture remains the same in a- and P-tubulin, still 
they differentiate themselves both sequential- 
ly and structurally. Sequentially they share 
40% residual identity, while structurally two 
positional gaps were observed in (3-tubulin at 
H1-S2 loop (45-46) and S loop (361-368). 
The larger gap in (3-tubulin gets well accom- 
modated by taxol Both the subunits get 
well associated with GTP; however, hydroly- 
sis is restricted to (3-subunit resulting in GTP 
formation at the exchangeable site. However, 
it gets sequestered in non-exchangeable site in 
a subunit ^^'^^ 

Tubulin with its active role in cell division 
has been considered as a potential anticancer 
drug target It comprises three distinct drug 
binding domains which include taxol, vinca- 
alkaloid and colchicine binding sites. Drugs 
associated with these sites were religiously 
involved in arresting the mitotic spindle for- 
mation Further, taxol and colchicine share 
overlapping residual interactions towards the 
inner surface of the microtubules Paclitaxel 
being the powerful drug for treating several 
solid tumors including breast, ovarian and 
non-small cell lung carcinomas prefers to 
bind with M-loop proximal to S loop resulting 



in stabilization of lateral interactions of two 
adjacent protofilaments The important res- 
idues associated with taxol binding include 
V23, D226, H229, T276, R278, R369 and 
Gly370. These residues were scattered around 
H1-S2 loop, H7, M-loop and S-loop ^\ Apart 
from paclitaxel, epothilones A and B, eleuthe- 
robin and discodermolide have also been re- 
ported to bind to the taxol binding site 

In spite of all these drug interactions, there 
is a report of drug resistance among these 
drugs. The main reason being cited is the re- 
sidual substitutions associated with the drug 
resistance. Even though literature supports the 
correlation between the residual substitutions 
and drug resistance, but still there is a contro- 
versy over their role because of the inclusion 
of pseudogenes during drug resistance analy- 
sis Thus, the debate on the role of residues 
on drug resistance remains elusive. The re- 
ported residual mutations in human P-tubulin 
include D26E, V60A, SI 72 A, PI 73 A, D197N, 
E198G, A231T, L240I, F270V, T274P, R282Q, 
Q292E, R306C, K350N, A364 and Y442C. 
Several types of resistance have been ob- 
served for the drugs like taxol, epothilone, 
hemiasterlin, 2-methoxyestradiol, vinca alka- 
loids and indancocine 

Previous studies report about drugs resist- 
ance and the reduced binding affinity of the 
available drugs, but still not much has been 
discussed about the drugs targeted against the 
mutant proteins. In this paper, we tried to tar- 
get the available chemical compounds from 
higher plants, lower plants and seaweed sec- 
ondary metabolites against wild and the mu- 
tants of human (3-tubulin to investigate their 
level of interactions. 

Materials and Methods 

To begin with, human (3-tubulin sequence 
was downloaded from SwissProt database 
(accession number: Q9BVA1.1) With no 
reported crystal structure of human (3-tubulin 
till date, a BLAST-PDB based search was 
carried out to identify suitable templates. Out 
of the reported hits, IJFF was downloaded 
from Protein Data Bank (www.rcsb.org) 



Auicenna Journal of Medical Biotechnology, \/oL d. No. 2, April-June 2014 



Selvaabumar C and Mohammed SMM 



which was further considered as a template 
(IJFF-B chain) for the modeUng of the query 
sequence. Pairwise ahgnment of the template 
and the query sequence were generated using 
Modeller9vl0 (Figure 1). Using single 
protocol template from Modeller, 30 struc- 
tures were generated. Of these generated 
structures, the model with the least DOPE 
score (Discrete Optimized Protein Energy) 
was considered for energy minimization with 
100 iterations using steepest descent in 
SwissPdb Viewer The template and the 
modelled structures were superimposed with 
each other (Figure 2A-C). The DOPE score of 
the template was -49714.55 and the modelled 
structure was -55870.25. The optimized struc- 
ture was further validated using PROCHECK 
of SAVES server (http://nihserver.mbi.ucla. 



edu/SAVES/) Similar methods were fol- 
lowed for the generation of mutants F270V, 
A364T and Q292E, respectively. 

Conversely, chemical compounds needed 
for docking against wild and the mutants of 
human (3-tubulin were obtained from marine 
flora. The chemical compounds were obtained 
from seaweed secondary metabolite database 
(www.swmd.co.in) All these available 
compounds were isolated form red algae Lau- 
rencia obtuse (RL) and Galaxaura marginata 
(RG type). Some of these compounds were 
reported to be cytotoxic against cancer cell 
lines During our previous study, out of 
the 517 compounds, we could identify four 
lead compounds named RG012 (6(3, 248- 
Dihydroxycholesta-4, 25-dien-3-one), RL381 
(diterpenes containing bromine), RL366 (15- 



Ijff B -REIVHIQAGQCGNQIGAKFWEVISDEHGIDPTGSYHGDSDLQLERINVYYNEAAGNKYVPRAILV 

humanbet a MREi vhiqagqcgnqigakfwevi sdehgi dptgtyhgdsdlqldri svyyneatggkyvprai lv 

Ijff B dlepgtmdsvrsgpfgqifrpdnfvfgqsgagnnwakghytegaelvdsvldwrkesescdclqg 

humanbet a dlepgtmdsvrsgpfgqifrpdnfvfgqsgagnnwakghytegaelvdsvldwrkeaescdclqg 

Ijff B fqlthslgggtgsgmgtlliskireeypdrimntfswpspkvsdtwepynatlsvhqlventde 

humanbet a FQLTHSLGGGTGSGMGTLLISKIREEYPDRIMNTFSWPSPKVSDTWEPYNATLSVHQLVENTDE 

Ijff B TYCIDNEALYDICFRTLKLTTPTYGDLNHLVSATMSGVTTCLRFPGQLNADLRKLAVNMVPFPRLH 

humanbet a TYCIDNEALYDICFRTLKLTTPTYGDLNHLVSATMSGVTTCLRFPGQLNADLRKLAVNMVPFPRLH 

1 j f f B F FM PG F AP LT S RG SQQY R ALT VP E LTQQM F DAK NMM AAC D P R HG R Y LT VAAVF RG RM SM K E VD EQM 

h um a n b et a F FM PG F AP LT S RG SQQY R ALT VP E LTQQVF DAK NMM AAC D P R HG R Y LT VAAVF RG RM SM K E VD EQM 

Ijff B LNVQNKNSSYFVEWIPNNVKTAVCDIPPRGLKMSATFIGNSTAIQELFKRISEQFTAMFRRKAFLH 

humanbet a LNVQNKNSSYFVEWIPNNVKTAVCDIPPRGLKMAVTFIGNSTAIQELFKRISEQFTAMFRRKAFLH 

1 j f f B WYTG EGM D EM E FT E AE S NM N D L VS E YQQYQD 

humanbet a WYTGEGMDEMEFTEAESNMNDLVSEYQQYQDATAEEEEDFGEEAEEEA 

Figure L Pairwise sequence alignment of template (IJFF-B chain) with the query sequence of human P-tubulin 




Figure 2. 3D structure of the template; A) the modelled structure of human P-tubulin; B) and its structural superimposition C) 
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epi-prostaglandin A2 diester), and RL376 
(diterpenes containing bromine) which exhibit 
better binding to wild type human (3-tubuUn 
protein All these selected compounds were 
taken into consideration for docking against 
the wild and mutant human (3-tubulin proteins 
in the current study. 

Next, through literature survey, we identi- 
fied chemical compounds both from higher 
and lower plants. To begin with, the com- 
pounds from higher plants include Berb- 
amine, Butulinic acid, Camptothecin, Cucur- 
bitacin, Ellipticine, Flavopiridol, Homohar- 
ringtonine, Silvestrol, Berberine, Daphnoretin 
and Podophyllotoxin. Regarding their sources, 
Berbamine is extracted from Berberis vulgar- 
is with a reported apoptosis in human myelo- 
ma cells Betulinic acid is a pentacyclic 
triterpenoid with reported antiretroviral, anti- 
malarial and anti-inflammatory activity and 
anticancer properties extracted from the bark 
of Betula pubescens 

Camptothecin isolated from the bark and 
stem of Camptotheca acuminate is a cytotoxic 
quinoline alkaloid which inhibits the DNA 
enzyme topoisomerase I Cucurbitacins are 
from the family of Cucurbitaceae with anti- 
cancer and anti-inflammatory activities 
Podophyllotoxin is a non-alkaloid toxin lig- 
nan extracted from the roots and rhizomes 
of Podophyllum species which is again an an- 
titumor agent Silvestrol is isolated from 
Aglaia foveolata and results in apoptosis in 
the cell lines of hormone-dependent human 
prostate cancer Homoharringtonine is ex- 
tracted from Cephalotaxus harringtonia 
which is identified to be a cytotoxic alkaloid 
and is generally reported to block the progres- 
sion of cells from Gl phase in to S phase and 
G2 phase into M phase Flavopiridol is an 
indigenous plant from India, which can arrest 
cell cycle progression at the Gl/S and G2/M 
boundaries Daphnoretin obtained from 
Wikstroemia indica exhibits strong antiviral 
and anti-tumor activities with a report of cell 
cycle arrest in the G2/M phase Ellipticine 
is isolated from Apocyanaceae plants with 
established antitumor and anti-HIV activities 



with their limited toxic side effects and their 
lack of hematological toxicity 

Regarding lower plants, Aclarubicin is pro- 
duced by Streptomyces galilaeus actively 
used in the treatment of cancer Dauno- 
rubicin and its derivative, doxorubicinare an- 
titumor anthracycline antibiotics are produced 
by Streptomyces peucetius . Blasticidin is a 
potent antifungal and cytotoxic peptidyl nu- 
cleoside antibiotic from Streptomyces griseo- 
chromogenes which plays a significant role in 
controlling prokaryotic and eukaryotic cell 
growth ^\ Chartreusin is a potent antitumor 
agent with a mixed polyketide-carbohydrate 
structure produced by Streptomyces chartreu- 
sis Neothramycin has been isolated from 
Streptomyces MC916-C4 which is a potent 
antitumor antibiotic of the pyrrolo (1, 4) ben- 
zodiazepine group Pirarubicin is an 
anthracycline drug which has a diversified 
antitumor activity All these chemical 
compounds were downloaded from Pubchem 
database. The chemical structures were ob- 
tained from chemical book (www.Chemical- 
book.com) (Table 1). The summary of the 
chemical compounds of lower and higher flo- 
ras along with marine derivates are tabulated 
in table 2. 

All these selected compounds were charged 
using gasteiger charges available in CHIME- 
RA software Energy was minimized using 
PRODRG and hydrogen atoms were added to 
them. Then these compounds were considered 
for docking using iGEMDOCK software 
The active sites of taxol (V23, D226, H229, 
T276, R278, R369 and Gly370) and epothi- 
lone (H227, A231, T274, R276, R282 and Q292) 

were selected for docking (Figure 3). Next, 
only the residual mutants with deleterious ef- 
fects were considered for homology model- 
ling. This was identified using amino acid 
substitution (AAS) tools like Polyphen2 (poly 
morphism phenotyping) PANTHER (pro- 
tein analysis through evolutionary relation- 
ship) and I-Mutant 2.0 ^\ The list of residu- 
al substitutions which bring in deleterious ef- 
fects are tabulated (Table 3) (for detailed re- 
port refer) Out of ten deleterious sites iden- 
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Table 1 . Chemical structures of higher and lower plants along with marine compounds 

RL381 Homoharringtonine 




Blasticidin 



I 



0 NH^ WH 



Chartreusin 



Daunorubicin 




Neothramycin 



HO 




Pirarubicin 




O OH O 



*'^ J^-^ QH 0 ^ 



RL376 




CH, 

q/ CH3 H ^' 


0 




RL366 





Camptothecin 
Cucurbitacin 





Berberine 




Daphnoretin 



OH 



Taxol 
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Table 2. Summary of chemical compounds of lower and higher floras along with marine derivatives 



Chemical compounds 


Molecular weight {g/mol) 


xlogP3 


H-bond donor 


H-bond acceptor 


Aclarubicin 


811.86 


3.8 


4 


16 


Blasticidin 


422.43 


-5.2 


6 


7 


Chartreusin 


640.58 


2.1 


5 


14 


Daunorubicin 


527.51 


1.8 


5 


14 


Neothramycin 


262.26 


0.2 


2 


5 


Pirarubicin 


627.63 


2.7 


5 


13 


RL376 


418.52 


3.1 


1 


7 


RL366 


695.81 


8.6 


2 


7 


RL381 


478.58 


2.6 


2 


9 


RG012 


414.62 


4.6 


2 


5 


Berbamine 


608.72 


6.1 


1 


8 


Butulinic acid 


456.70 


8.94 


2 


3 


Camptothecin 


348.35 


1 


1 


5 


Cucurbitacin 


514.65 


2.7 


4 


7 


EUipticine 


246.30 


4.8 


1 


1 


Flavopiridol 


401.84 


3.3 


3 


6 


Homoharringtonine 


545.62 


0.8 


2 


10 


Silvestrol 


654.65 


1.6 


4 


13 


Berberine 


336.36 


3.6 


0 


4 


Daphnoretin 


352.29 


3.3 


1 


7 


Taxol 


853.90 


2.5 


4 


14 


Epothilone 


493.65 


4 


2 


8 


Podop hyllotoxin 


414.40 


2 


1 


8 



Table 3. Predicting the effect of residual substitutions through amino acid substitutions tools 
Polyphen2, PANTHER and I-Mutant2.0 



Residual position 




Polyphen2 


PANTHER 


I-Mutant 2.0 




Prediction score 


Pdeleterious 


(in kcal/mot) 


V60A 




0.88 


Possibly damaging 


0.81 


-2.2 


P173A 




NA 


NA 


0.67 


-1.9 


D197N 




0.99 


Probably damaging 


0.84 


-0.37 


L240I 




0.95 


Probably damaging 


0.68 


-1.12 


F270V 




0.9 


Possibly damaging 


0.5 


-1.66 


Q292E 




0.99 


Probably damaging 


0.6 


0.28 


R306C 




1 


Probably damaging 


0.9 


0.1 


K350N 




0.99 


Probably damaging 


0.79 


-0.75 


A364T 




0.88 


Probably damaging 


0.76 


-1.3 


Y422C 




1 


Probably damaging 


0.92 


1.54 



tified, only three were selected as these sub- 
stitutions were found to be proximal to the 
drug binding sites. Unfortunately, the rest of 
the positions were located at distal sites of 
both taxol and epothilone binding sites. Dock- 
ing of the lower and higher plants along with 
marine isolates were carried out using iGEM 



DOCK for taxol-wild, taxol-mutants (F270V, 
A364T), epothilone-wild and epothilone- 
mutant (Q292E). The docking study followed 
the accurate docking protocol which was very 
slow with population size of 800 and the 
number of solution equal to 10. The genera- 
tion number was maintained at 80. Similar 
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Figure 3. Drug binding sites in human |3-tubulin for epo- 
thilone (Red) and taxol (Blue) 

protocol was followed for taxol and epothi- 
lone binding pockets. Finally, we measured 
the size of the normal grooves for wild and 
the mutants of taxol and epothilone binding 
sites using SwissPdb Viewer. 

Results 

In the present study, the chemical com- 
pounds obtained from different resources like 
lower plants, higher plants and marine isolates 
were docked against the wild and the mutant 



human (3-tubulin proteins. This study was 
mainly instigated to identify the level of inter- 
actions exhibited by both the wild and the 
mutants for taxol resistance (F270V, A3 64) 
and epothilone resistance (Q292). For this 
study, we considered the results of our previ- 
ously reported lead molecules from seaweed 
secondary metabolites. Simultaneously, the 
chemical compounds from lower and higher 
plants were also taken into consideration. 

First of all, the wild tubulin when docked 
with these compounds exhibited better bind- 
ing to Aclarubicin. Next, the compounds from 
higher plants like Berbamine, Butulinic acid, 
Cucurbitacin, taxol and Podophyllotoxin ex- 
hibited better interactions in the wild type in 
the taxol binding site similar to the marine 
compounds. Next, the mutants were consid- 
ered for docking with all the available com- 
pounds. Here, Campothecin, Flavopiridol and 
Berberine exhibited better interactions in mu- 
tantl (F27 OV). This site is proximal to M- 



Table 4. Docking of chemical compounds from lower, higher and marine resources against taxol and epothilone binding sites in 

wild and mutant proteins 



Origin 



Compounds 



Taxol binding site 
(binding energy in kcal/mol) 



Epothilone binding site 
(binding energy in kcal/mol) 





Wild 


Mutantl (F270V) 


Mutant! (A364T) 


Wild 


Mutant (Q292E) 


Aclarubicin 


-125.84 


-121.20 


-96.84 


-69.66 


-69.68 


Blasticidin 


65.17 


13.35 


88.07 


-33.98 


72.76 


chartreusin 


319.11 


38.04 


388.28 


221.045 


54.33 


Daunombicin 


207.93 


-4.9 


46.49 


0.45 


-11.00 


Neothramycin 


-66.27 


-66.98 


69.74 


-28.26 


-31.36 


Pirambicin 


90.65 


95.46 


50.78 


332.44 


66.66 


RL376 


-86.97 


-83.51 


-81.81 


-45.38 


-33.72 


RL366 


-82.61 


-79.19 


-81.44 


-30.86 


-37.99 


RL381 


-78.16 


-75.38 


-79.92 


-28.33 


-30.66 


RG012 


-77.55 


-76.07 


-65.17 


-38.40 


-38.49 


Berbamine 


-19.55 


235.23 


129.56 


26.44 


-33.49 


Butulinic acid 


-79.00 


-17.77 


-15.62 


-20.91 


2.02 


Camptothecin 


-82.97 


-79.27 


-86.75 


-44.01 


-43.85 


Cucurbitacin 


-86.92 


-84.52 


-78.20 


-24.33 


-24.36 


Ellipticine 


-66.21 


-64.62 


-75.40 


-39.95 


-15.65 


Flavopiridol 


-67.46 


-69.65 


-74.93 


-4.32 


-39.29 


Homoharringtonine 


107.36 


145.86 


13.95 


181.80 


7.39 


Silvestrol 


133.96 


153.26 


57.28 


43.01 


25.84 


Berberine 


-75.72 


-81.20 


-82.74 


-29.36 


-36.08 


Daphnoretin 


-80.37 


-79.21 


-82.75 


-34.01 


-37.03 


Taxol 


-82.60 


-79.20 


-81.48 






Epothilone 








-37.10 


-21.84 


Podophyllotoxin 


-83.00 


-80.23 


-75.71 


-40.19 


-38.18 



Lower plants 



Marine 
compounds 



Higher plants 
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loop and is also the taxol binding site. With 
respect to mutant2 (A364T), the chemical 
compounds like Camptothecin, Ellipticine, 
Flavopiridol, Berberine, Daphnoretin showed 
better chemical interactions. In epothilone 
binding site, lower plant compounds like 
Aclarubicin and Blasticidin showed better in- 
teraction in the wild type. The compounds 
like Neothramycin, Daunorubicin and Aclaru- 
bicin showed slightly better interaction with 
the mutation of glutamate with Glu at residual 
position 292. Thus, our overall study confirms 
that marine compounds exhibit better interac- 
tion in this mutant protein except RL376. 
Among higher plant products, Camptothecin, 
Ellipticine, Epothilone, Podophyllotoxin ex- 
hibit better interaction in the wild type of 
epothilone binding site. But compounds like 
Berbamine, Flavopiridol, Berberine and Daph- 
noretin conveyed better interactions again in 
the mutant (Table 4). General survey confirms 
a similar type of interactions exhibited both 
by wild and the mutant proteins. Thus, a clos- 
er inspection of the pocket was much needed 
for both sets of proteins. Therefore, we meas- 
ured the groove size for wild, F270V and 
A364T of the taxol binding site. The volume 
was 550, 607 and 220 respectively. How- 
ever, for epothilone binding pockets, the wild 
and the mutant (Q292E) maintained same 
volume of 8 A^. There was an increase in 
pocket volume observed for mutant 1 in com- 

A 



parison with the wild type due to non- 
conservative residual substitution. Similarly, 
the volume of mutant2 decreases in compari- 
son with the wild type in spite of conservative 
substitution. However, in epothilone binding 
site, the conservative substitution does not 
bring in any change in the volume both in 
wild and the mutant and the grooves were 
found to be discontinuous with the maximum 
volume size of 8 A^ (Figure 4). 

Discussion 

During the docking of the wild and the mu- 
tants against the available compounds, only 
three mutants were generated which were 
proximal to the drug binding site, and were 
also reported to be lethal. Thus, the generated 
mutants were christened as mutant 1 (F270V) 
mutant2 (A364T) for the taxol binding site. 
Similarly, mutant (Q292E) was selected for 
epothilone binding site. Docked structures 
were separately investigated for their interac- 
tions using CHIMERA software. We observed 
that Neothramycin exhibits interactions with 
the taxol binding site through P272 and R276 
(Figure 5A). Similar interactions were exhib- 
ited by mutant 1 (Figure 6A). Also mutant2 
exhibited interaction with L361 (Figure 6B). 
Likewise, RL381 displayed interactions in the 
wild type of taxol binding site with H227, 
R282, and R359 (Figure 5B). Camptothecin 
shows interaction with T274 (Figure 5C) 




5 4T^''1-< 





Figure 4. The groove of the active site pockets in wild and the mutant human |3-tubuHn; A) wild-taxol binding site; B) mutantl -taxol 
binding site; C) mutant2-taxol binding site; D) wild-epothilone binding site; E) mutant-epothilone binding site 
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Figure 5. The docking of chemical compounds in the taxol binding sites of wild human p-tubulin; A) Neothramycin-taxol-wild; B) 
RL381-taxol-wild; C) Camptothecin-taxol-wild 



while taxol based mutants, irrespective of 
their higher binding energy, do not display 
any interactions. Ellipticine also confirms no 
significant binding in mutant2. Flavopiridol 
does not show any significant interactions in 
the wild taxol binding site. But mutant 1 
shows interaction with H227 (Figure 6C). Al- 
so mutant2 shows interactions with T274 and 
A275 (Figure 6D). Berberine also shows no 
significant interaction in the wild type and 
mutant2 of taxol binding site but mutant 1 
confirms their interaction with R276 (Figure 
6E). Again with Daphnoretin, both wild and 
mutant2 show no interactions, but mutant 1 
interacts with R276 (Figure 6F). Next, re- 
garding the epothilone binding site, Neothra- 
mycin binds with mutant at T274 (Figure 7A). 
RL366 interacts in the mutant with R276 



(Figure 7B). However, the rest of the com- 
pounds like RL381, Flavopiridol, Berberine 
and Daphnoretin exhibited no significant in- 
teractions. Since no remarkable change in the 
pattern of interactions both in wild and the 
mutant was observed, we considered the 
groove analysis for the taxol and the epothi- 
lone binding site. This study confirms an in- 
crease in the channel size for mutant 1 and a 
decrease in the channel size for mutant2 in 
comparison with wild type. Moreover, epothi- 
lone binding site remains undisturbed with the 
residual substitution. We hypothesized that 
the drug resistance especially for taxol bind- 
ing drugs could be related to the conforma- 
tional changes in the active site pocket due to 
residual mutations. 




Figure 6. The docking of chemical compounds in the taxol binding sites of mutant human P-tubulin. A) Neothramycin-taxol-mutantl; 
B) Flavopiridol-taxol-mutantl; C) Berberine-taxol-mutantl; D) Daphnoretin- taxol-mutantl; E) Neothramycin-taxol-mutant2; 
F) Flavopiridol-taxol-mutant2 
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Figure 7. The docking of chemical compounds in the epothilone binding sites of mutant (Q292E) human P-tubuhn. A) Neothramycin- 
epothilone-mutant; B) RL366-epothilone-mutant 



Conclusion 

The residual interaction analysis in wild 
and the mutants of taxol and epothilone bind- 
ing sites reveals a better drug binding with the 
lower plants in the wild types of taxol and 
epothilone. However, the marine compound 
and the higher plants in spite of their increase 
in binding score could not establish a better 
contact with the residues of both the binding 
sites. Further, the chemical compounds show- 
ed interactions with the charged residues like 
arginine and histidine for better drug binding. 
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